Magnetic moments of vector, axial, and tensor mesons in lattice QCD 
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We present a calculation of magnetic moments for selected spin-1 mesons using the techniques 
of lattice QCD. This is carried out by introducing progressively small static magnetic field on the 
lattice and measuring the linear response of a hadron's mass shift. The calculations are done on 24* 
quenched lattices using standard Wilson actions, with /3=6.0 and pion mass down to 500 MeV. The 
results are compared to those from the form factor method. 
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I. INTRODUCTION 

Magnetic moment is a fundamental property of 
hadrons that arises from the linear response of a bound 
system to an external stimulus (in this case, magnetic 
field). It is a good testing ground for studying the in- 
ternal structure of hadrons as governed by the quark- 
gluon dynamics of QCD, the fundamental theory of the 
strong interaction. Efforts to compute the magnetic 
moment on the lattice come in two categories. One 
is the form factor method which involves three-point 
functions P, i, i, i, H, i, 0- The other is the back- 
ground field method using only two-point functions (mass 
shifts) is, ^, J^, Jd]. The form factor method requires an 
extrapolation to zero momentum transfer Gm{Q^ = 0) 
due to the non-vanishing minimum discrete momentum 
on the lattice ^M,- The background field method, on 
the other hand, accesses the magnetic moment directly 
but is hmited to static properties due to the use of a 
static field. Here we report a calculation of the vector 
meson magnetic moments in this method, in parallel to 
a recent calculation in the form factor method [13,]. We 
also report results mesons in the axial and tensor sectors. 
Some of the preliminary results have been reported in a 
conference |14| . It is an extension of our earlier work 
on baryon magnetic moments [1^1 and electric [la| and 
magnetic polarizabilities [l7| in the same method. 



II. CORRELATION FUNCTIONS 

The mass of a meson can be extracted from the time- 
ordered, two-point correlation function in the QCD vac- 
uum, projected to zero momentum 



G{t)^Y.^^{xW{Q)) 



(1) 



where rj is the interpolating field of the meson under con- 
sideration. The general form of the interpolating field 
with a simple qiq2 quark content can be written as 



where F is a general gamma matrix that depends on the 
meson type. For the spin-1 mesons considered in this 
work, F = 7p for vector mesons, F = 757/j for axial 
mesons, and F = 7^7j/ for tensor mesons. Here we only 
consider mesons with gi and (j2 different to avoid the 
complication of disconnected loops. 

On the quark level, Eq. ([T]) is evaluated by contracting 
out the quark pairs, 

G{t) = - 5] Tr [S,, {x, 0horhol5Sl {x, Oh^T] (3) 

X 

where ^^(a:, 0) denotes the fully- interacting quark prop- 
agator. It is defined as Euclidean-space path integrals 
over gauge field G^ 



Sgix,0) 



J DG^det{M)e~'^^M- 
J DG^det{M)e-So 



(4) 



where Sq is the gauge action of QCD and M = ^^D^ -\- 
niq its quark matrix. On the lattice, the propagator is 
evaluated numerically by Monte-Carlo methods. We use 
quenched approximation in this work which corresponds 
to setting det(M ) to a constant. 

On the hadronic level, the correlation function is sat- 
urated by the complete spectrum of intermediate states 
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(5) 



where mi are the masses and Wi are spectral weights that 
are a measure of the ability of the interpolating field to 
excite or annhilate the states from the QCD vacuum. 
The ground state can be extracted by fitting G(t) at large 
time. 

To compute magnetic moments, we need to use polar- 
ized interpolating fields. For a magnetic field applied in 
the z-direction, we use 



v± = —^qi {t'^x 



i^y)'l2 = —^iVx±iVy)- (6) 



V = qiTq2 



(2) 



The interaction energies E± are extracted from the cor- 
relation functions 

('7±'7±> = 2 [^'^^''^^ ^ * {{Vxvl) - iVyvD) + iVyVl)] ■ 

(7) 



Eq. ([7]) implies that the polarization comes from the 
imaginary parts of the off-diagonal correlation between x 
and y components in the presence of the magnetic field. 
These imaginary parts are zero in the absence of the field, 
so they are responsible for the magnetic moments we ob- 
serve. We use for vector mesons Tx = 71 and Fj, = 72; 
for axial mesons T^ = 7571 and Ty = 7572; and for tensor 
mesons T^: = 7273 and Ty = 7173. 

For each meson type, different quark combinations gi 
and 92 correspond to different states. In the case of vector 
mesons, they are the well-known p'^(du), p"(ud), (jjiss), 
K*+(su), K*~{us, and K*°{sd). In the case of axial 
mesons, they are ai{du), a^{ud), Kl'^{su), K^~{us), 
and Kl^{sd). In the case of tensor mesons, they are less 
well-known and we call them bj{du), bi{ud), K^ {su), 
Kl~{us, and K'^^{sd). Counting the states with ss con- 
tent (like the meson), we cover 18 states of spin-1 
mesons. 



III. BACKGROUND-FIELD METHOD 



For a particle of spin s in uniform fields, 



E± = m ± fiB 



(8) 



where the upper sign means spin up and the lower sign 
means spin-down relative to the magnetic field, and /x = 
g-^s. We use the following method to extract the g 
factors. 



g = m- 



{E, 



{E- 



eBs 



(9) 



In order to place a magnetic field on the lattice, we con- 
struct an analogy to the continuum case. The covariant 
derivative of QCD is modified by the minimal coupling 
prescription 



D^ = d^ 



gG^ + qA^ 



(10) 



where q is the charge of the fermion field and A^ is 
the four- vector potential describing the background field. 
On the lattice, the gluon fields G^ are introduced via 
link variables U^{x) = eyi\i {igaG ^) . So the prescription 
amounts to multiplying a U(l) phase factor exp(iga^^) 
to the gauge links. Choosing Ay = Bx, a constant mag- 
netic field B can be introduced in the z-direction. Then 
the phase factor is applied to the y-links 



a new dynamical ensemble is needed that couples to u- 
quark (q=l/3), d-and s-quark (q=-2/3). This requires a 
Monte Carlo algorithm that can treat the three fiavors 
distinctively. Quark propagators are then computed on 
the ensembles with matching field values. This has not 
been attempted. The second can be termed as the re- 
weighting method in which a perturbative expansion of 
the action in terms of the field is performed. There has 
been an attempt [18[ to compute the neutron electric po- 
larizability in this method. It involves the evaluation of 
disconnected diagrams. The third is what we call U(l) 
quenched. No field is applied in the Monte-Carlo genera- 
tion of the gauge fields, only in the valence quark propa- 
gation in the given gauge background. In this case, any 
gauge ensemble can be used to compute valence quark 
propagators. 

We use standard Wilson actions on the 24'* lattice 
at /3 = 6.0, both SU(3) and U(l) quenched, and six 
kappa values k=0.1515, 0.1525, 0.1535, 0.1540, 0.1545, 
0.1555, corresponding to pion mass of about 1015, 
908, 794, 732, 667, 522 MeV. The critical value of 
kappa is Kc =0.1 571. The strange pion mass is set at 
K=0.1535. The source location for the quark propaga- 
tors is (x,y,z,t)=(12,l,l,2). We analyzed 100 configu- 
rations. The following five dimensionless numbers 77 = 
qB 0^^+0.00036, -0.00072, +0.00144, -0.00288, -h0.00576 
give four small B fields (two positive, two negative) at 
eBa^=-0.00108, -(-0.00216, -0.00432, -1-0.00864 for both 
u and d (or s) quarks. These field values do not obey 
the quantization condition for periodicity since the val- 
ues given by the condition cause too strong (too large a 
mass shift) for the small-field-expansion method to work. 
To minimize the boundary effects, we work with Dirich- 
let boundary conditions in the x-direction and large N^. 
In addition, we place the source in the middle of the lat- 
tice in the x-direction so that quarks have little chance 
of propagating to the edge. We also use Dirichlet bound- 
ary conditions in the t-direction to maximize the num- 
ber of time slices for mass extraction. To eliminate the 
contamination from the even-power terms, we calculate 
mass shifts both in the field B and its reverse —B for 
each value of B, then take the difference and divide by 
2. Another benefit of repeating the calculation with the 
field reversed is that by taking the average of dm{B) and 
6m{—B) in the same dataset, one can eliminate the odd- 
powered terms in the mass shift. The coefficient of the 
leading quadratic term is directly related to the magnetic 
polarizability |17| . 



U„ 



exp [iqaBx)Uy. 



(11) 



In our calculations, we use a linearized version for small 
field strengths 



Uy -^ {\ + i qaBx)Uy. 



(12) 



The computational demand of such background-field cal- 
culations can be divided into three categories. The first is 
a fully- dynamical calculation. For each value of the field. 



IV. RESULTS AND DISCUSSION 

A. Vector mesons 

Fig. [B displays a typical effective mass plot for /9+. 
Both the mass and the mass shifts are shown. Good 
plateaus exist for all six pion masses. The mass shifts 
are extracted from the time window 10 to 13, as indi- 
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FIG. 1: Effective mass plot for the p"*" vector meson mass 
at zero field (top), and effective mass shifts at the weakest 
magnetic field (bottom) in lattice units. The solid and empty 
symbols correspond to the heaviest and lightest pion masses, 
respectively. 



cated in the figure. Fig. [2] shows the mass shifts, defined 
a.s S — g{eBs) from Eq. ([5]), as a function of the field 
for the p+ meson. The slope gives the g-factor. There 
is good linear behavior going through the origin at all 
the field values, an indication that contamination from 
the higher-power terms has been effectively eliminated by 
the {5{B) — S{—B))/2 procedure. This is also confirmed 
numerically by the smallncss of intercept as shown in the 
fit results y — ax + b. At the lightest pion mass, there 
is a slight deviation from linear behavior at the stronger 
fields. For this reason, we only use the two smallest field 
values to do the linear fit at all the pion masses. 

Fig. [3] shows the g-factors for the vector mesons as a 
function of pion mass squared. The lines are simple chiral 
fits using the ansatzs 



g = ao +aim^, 



and 



ao + aim„ 



02771^. 



(13) 



(14) 



probably desirable. But overall the g-factors have a fairly 
weak pion mass dependence. At large pion masses, the 
g-factor of p+ approaches 2, consistent with a previous 
lattice calculation using the charge-overlap method 13|. 
Our results for p+ are slightly higher than those from the 
form factor method (see Fig. 8 in Ref. [13]). The results 
confirmed that gp- = —gp+ and gx-- = ~9k*+- We also 
confirmed gpO = numerically (not shown). These re- 
lations are expected from symmetries in the correlation 
functions (these particles are charge eigenstates) . The re- 
sults also show that as far as g-factors are concerned the 
p mesons arc quite similar to their strange counterparts 
K* mesons. 

Note that the extracted g-factors are in the particle's 
natural magnetons. To convert them into magnetic mo- 
ments in terms of the commonly-used nuclear magnetons 
(pat), we need to scale the results by the factor 938/M 
where M is the mass of the particle measured in the same 
calculation at each pion mass. Fig. |4] shows the results 
for p+ and K*^. The different pion-mass dependence be- 
tween p+ and K*^ mostly comes from that in the masses 
that are used to convert the g-factors to magnetic mo- 
ments. The values at the chiral limit extrapolated from 
Eq. do]) are ^p+ = 3.25(3)/XAr and ij.k*+ = 2.81(1)/X7v. 
There is no experimental information on these quantities. 
Compared to the form factor method (see Fig. 7 in T3]), 
our results are again a little higher. At the strange pion 
mass point (the 3rd data point from the left), the two 
coincide to give a prediction for the magnetic moment of 
the (/)(1020) meson, p^ = 2.07(7)pAr. 

Fig. m shows the results for K*'^. Our results confirm 
the expectation that fJ.K*o is small but has an interesting 
quark mass dependence. It is positive when the d-quark 
is heavier than the s-quark, exactly zero when they are 
equal, and turns negative when the d-quark is lighter 
than the s-quark. The same behavior has been observed 
in the form factor method (see Fig. 11 in [13|). 



They serve to show that there is onset of non-analytic be- 
havior as pion mass is lowered, so a linear extrapolation is 
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FIG. 2: Mass shifts for the p^ meson as a function of the magnetic field in lattice units at the six pion masses (heavy to 
light from top left to right, then to bottom left to right). The slope of the mass shift at each pion mass gives the g factor 
corresponding to that pion mass. The line is a fit using only the two smallest B values. 
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FIG. 3: G-factors for the p (top) and K* (bottom) vector 
mesons as a function of pion mass squared. The 2 hues are 
chiral fits according to Eq. (fT3)) (dashed), Eq. ([T4|l (dotted). 
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FIG. 5: Magnetic moments (in nuclear magnetons) for K* . 
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4: Magnetic moments (in nuclear magnetons) for p^ and 



B. Axial mesons 



Fig. [HI shows the effective mass shifts for the af axial 
meson at the 2nd value of the magnetic field {eBa^ = 
0.00216). The signal is noisier compared to the vector 
case, but a plateau is still visible between time slice 3 to 
5 in the mass shifts. Fig.[7]shows the g- factors for af and 
K^ extracted from this window at the six pion masses. 
The g-factors are very similar to their counterparts in the 
vector channel (see Fig. ^. Fig. [5] shows the g-factors 
for K^'^. They are small as expected, but have a linear 
behavior across the zero as a function of the pion mass 
squared. Interestingly, they have the opposite sign to 
that in the vector channel (see Fig. [5]): negative when the 
d-quark is heavier than the s-quark, exactly zero when 
they are equal, and turns positive when the d-quark is 
lighter than the s-quark. 
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FIG. 7: G-factors for the af (top) and K^"^ (bottom) axial 
mesons as a function of pion mass squared. The 2 lines are 
chiral fits according to Eq. (fT5)) (dashed), Eq. ^^ (dotted). 
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FIG. 8: G-factor for the neutral axial meson Ki . 



FIG. 6: Effective mass plot for the af axial meson mass at 
zero field (top), and effective mass shifts at the 2nd weak- 
est magnetic field (bottom) in lattice units. The solid and 
empty symbols correspond to the heaviest and 2nd lightest 
pion masses, respectively. 



Tensor mesons 



Fig. [9] shows the effective mass shifts for the bf tensor 
meson at the 2nd value of the magnetic field [eBa? = 
0.00216). There is a signal, but much noisier than the 
axial case and there is barely a plateau in the mass shifts. 



If we fit the data between time slice 4 to 6, the results 
are shown in Fig. [TO] for bf and K^ extracted from this 
window at the six pion masses. The g-factors have a weak 
pion mass dependence and relatively large errors. They 
have smaller values than the axial counterparts. Fig. [TT] 
shows the g-factors for the neutral Kf'^, which displays a 
similar linear behavior as the axial counterpart. This is 
evidence that there is indeed a signal in the tensor case. 
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FIG. 10: G-factors for the charged tensor mesons h^ (top) 
and Kl (bottom) as a function of pion mass squared. The 
2 lines are chiral fits according to Eq. (|13p (dashed) , Eq. (|14p 
(dotted). 



FIG. 9: Effective mass plot for the \y[ tensor meson mass at 
zero field (top), and effective mass shifts at the 2nd weak- 
est magnetic field (bottom) in lattice units. The solid and 
empty symbols correspond to the heaviest and 2nd lightest 
pion masses, respectively. 




FIG. 11: G-factor for the neutral tensor meson K* 



V. CONCLUSION 

In conclusion, wc have computed the magnetic mo- 
ment of vector, axial and tensor mesons on the lattice 
using the background field method and standard lattice 
technology. Our results for the vector mesons are con- 
sistent with those from the form factor method where a 
comparison is possible. The results for the axial and ten- 
sor mesons are new, although the latter still suffer from 
large errors. Nonetheless, our results demonstrate that 
the method is robust and relatively inexpensive. Only 
mass shifts are required. There is no experimental in- 
formation on these quantities so the lattice results can 
serve as a guide from first principles. Since the feasibil- 
ity of the method is extended to the meson sector, the 
calculation can be improved in a number of ways. First, 
it should be repeated on a lattice of larger size in order 
to get an idea about finite-volume effects. More statis- 
tics are needed in the axial and tensor cases to better 
isolate the signals. Second, there is a need to push the 



calculations to smaller pion masses so that reliable chi- 
ral extrapolations can be applied. Third, the calculation 
should be extended to full QCD in order to see the effects 
of the quenched approximation, both in the SU(3) sector 
and in the U(l) sector. With the availability of dynami- 
cal configurations, all of the improvements can be made 
at the same time. In particular, the U(l) effect in the 
sea quarks can be evaluated by re-weighting the deter- 
minants in the correlation functions, without the need to 
generate new dynamical ensembles. 
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